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Abstract 

In this paper, the problem of joint design of Spectrum Sensing (SS) and receive beamforming 
(BF), with reference to a Cognitive Radio (CR) system, is considered. The aim of the proposed design 
is the maximization of the achievable average uplink rate of a Secondary User (SU), subject to an 
outage-based Quality-of-Service (QoS) constraint for primary communication. A hybrid CR system 
approach is studied, according to which, the system either operates as an interweave (i.e., opportunistic) 
or as an underlay (i.e., spectrum sharing) CR system, based on SS results. A realistic Channel State 
Information (CSI) framework is assumed, according to which, the direct channel links are known by the 
multiple antenna receivers (RXs), while, merely statistical (covariance) information is available for the 
interference links. A new, closed form approximation is derived for the outage probability of primary 
communication, and the problem of rate-optimal selection of SS parameters and receive beamformers 
is addressed for hybrid, interweave and underlay CR systems. It is proven that our proposed system 
design outperforms both underlay and interweave CR systems for a range of system scenarios. 
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I. Introduction 

Spectrum scarcity, as it had been observed in 2002 by the Federal Communications Com¬ 
mission (FCC) BH, constitutes a major drawback, in terms of facilitating wireless communi¬ 
cations services. To overcome such an obstacle, the notion of Cognitive Radio (CR) ll^-[[6l 
was introduced, targeting at improving the information throughput by optimally exploiting the 
under-utilized spectrum. 

In practice, two different categories of CR systems have been devised: a) Underlay (or 
spectrum sharing) CR systems, where a Primary User (PU) allows the reuse of its spectrum 
by an unlicensed Secondary User (SU), provided that the interference received by the PU will 
be such that an interference temperature constraint will not be violated, and b) Interweave (or 
opportunistic) CR systems, where the SU senses the spectrum environment and transmits at 
time intervals during which primary activity is not detected. As it has been explained in [[71, 
via an analytical comparative study, each of the described CR approaches is characterized by 
drawbacks of different kind. For instance, the throughput performance of an interweave system is 
seriously affected by the quality of Spectrum Sensing (SS), while underlay CR systems, in turn, 
manipulate their transmission strategy according to a fixed interference temperature constraint, 
without exploiting the traffic pattern (or activity profile) of the PU. 

With the aim of relaxing such inherent drawbacks, a hybrid interweave/underlay CR approach 
has been investigated in the literature, in order to exploit the benefits of the two standard CR 
approaches. However, the full potential of such a scheme, considering a realistic and practical 
system, from a Channel State Information (CSI) viewpoint, has not been studied so far to the 
best of our knowledge. For instance, in works such as [fSll- lfm . hybrid CR systems are proposed, 
however, either no average rate-based performance analysis under channel fading is undertaken 
[|8l , [[Toll or the unrealistic assumption of perfect SS is assumed BH. In BITI a Single-Input- 
Single-Output (SISO) framework is investigated, thus, not being in accordance with today’s most 
wireless systems, where multiple antennas are used at the Base Stations (BSs) and (possibly) 
at the mobile devices as well. Furthermore, in BT2ll . the problem of joint, optimal (in terms 
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of average SU rate) SS and power poliey design is investigated for a hybrid CR system in 
the uplink, however, assuming the existenee of uneorrelated reeeive antennas and applying a 
Maximal Ratio Combining (MRC) reeeiver. Also, in lfT3l . the downlink of a Multiple-Input- 
Multiple-Output (MIMO) hybrid CR system is studied analytieally and performanee eomparisons 
are made with the standard interweave and underlay CR systems. Nonetheless, the existenee of 
spatially uneorrelated antennas at the transmitters (TXs) is assumed, along with the applieation 
of a simple, truneated power alloeation seheme, depending on an interferenee temperature 
threshold. Moreover, in m, the problem of optimal, in terms of the aehievable average uplink 
rate, beamforming (BF) problem is presented and solved, foeusing on the two-user, multiple- 
antenna interferenee ehannel, with eombined instantaneous and statistieal CSI. However, the two 
systems are eharaeterized by the same priority, thus, no solution for the equivalent CR system 
was provided. 

Motivated by the above, in this paper we foeus on the uplink of a hybrid interweave/underlay 
CR system. The hybrid CR system operates either as an interweave or as an underlay CR system, 
based on the results of the SS proeedure. In sueh a setting, our eontributions ean be summarized 
as follows: 

• Foeusing on a spatially eorrelated fading ehannel model and assuming a eombined CSI 
setting at the reeeivers (RXs) (CSIR), where direet links are known instantaneously and 
interferenee links are merely known based on their seeond order statisties, we derive 
new elosed form approximations for the outage probability of primary eommunieation, 
eonsidering the hybrid CR system as well as the standard interweave and underlay CR 
systems. Simulations show that the derived expressions approximate the aetual outage 
probability suffieiently well. 

- Foeusing on primary systems applying MRC reeeivers, the derived approximations are, 
to the best of our knowledge, the first appearing in the literature, that, unlike works sueh 
as llT4l . also inelude additive noise and do not presume a speeifie relation between the 
eovarianee matriees of the desired and interfering ehannels. Also, in eontrast with [fT5l . 
both the desired and the interferenee links are spatially eorrelated. 

• Having derived the expressions deseribed above, and foeusing on an interferenee-limited 
system, i.e., a system for whieh interferenee is the dominant souree of signal degradation. 
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as compared to noise, for the first time, we formulate and solve the problem of jointly 
determining a) the transmit power of the SU, b) the applied receive BF scheme, as well 
as c) the SS parameters, such as to maximize the achievable average rate of the SU, 
subject to an outage-based constraint on primary communication. The derived optimization 
framework is applied to all previously described CR system approaches, i.e., hybrid, as well 
as interweave/underlay. 

- The derived optimization framework can be applied for determining the transmit power 
and optimizing the BF and SS design for uplink communication of CR systems as well as 
for Licensed Shared Access (LSA) systems ifT^ . [fTTIl . where the operation of a licensee 
user without violating the performance of an incumbent user, is crucial. 

- Focusing, in particular, on the SS and receive BF optimization framework, we note that, 
to the best of our knowledge, the BF and SS problems are treated in a joint manner, for 
the first time. 

• The throughput performance of the optimized hybrid CR system is evaluated and compared 
to the performance achieved by the two optimized standard CR systems. It turns out that 
the hybrid system outperforms the standard ones for the whole range of values of the 
investigated system design parameters, i.e., the outage constraint and the activity profile of 
the PU. It is also shown that the performance of the hybrid CR system for low primary 
activity profiles, converges to the one achieved by the interweave system, while, for high 
primary activity profiles, the hybrid CR system behaves in a similar manner as the standard 
underlay one. 

The following notations are adopted throughout the paper: all lower case boldface letters 
indicate vectors, whereas all upper case boldface letters denote matrices. Superscript stands 
for Hermitian transpose, || ■ || denotes the Euclidean norm and Pr{A) denotes the probability 
of event A. Symbol denotes the {p,q)-th element of matrix A. The all-zero vector of 

dimension nx 1 is denoted as 0„. The identity matrix of dimension nxn is denoted as !„, whereas 
E|x{/(A, y)} symbolizes the conditional (with respect to Random Variable (RV) X) expectation 
of function f{X, Y). Also, tr(A), Aj(A) and rank(A) denote the trace, the j-th largest eigenvalue 
of square matrix A and its rank, respectively. For a random vector x,x ^ CM {pi ^ S) denotes 
that X follows a Circularly Symmetric Complex Gaussian (CSCG) distribution, with mean /x 
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and covariance matrix S. Furthermore, exp(-) and In(-) denote the exponential and logarithmie 
funetions. Additionally, -Ei(-) represents the exponential integral funetion, as defined in [fT^ eq. 
(5.1.1)] and Q(-) represents the eomplementary Gaussian distribution funetion, as defined in [fT9l 
eq (4.1)]. Finally, 7 ~ 0.5772 stands for the Euler-Maseheroni eonstant, as defined in [|T8l eq. 
(4.1.32)]. 

II. System Model 

A. Signal and channel model 

The uplink of a CR system is eonsidered, as shown in Fig. [H whieh eomprises of a single¬ 
antenna TX of a primary network, TX p, that eommunieates with a multiple-antenna RX, RX p. It 
is assumed that the primary network is willing to share part of its speetral resourees with a 
seeondary network. The latter is eomposed of a single-antenna TX, TX s, eommunieating with 
a multiple antenna RX, RX s. In what follows, it is assumed that RX p and RX s are equipped 
with M antennas, eaeh. 

The Single-Input-Multiple-Output (SIMO) ehannel between TX m and RX n is denoted as 
hmn e m,n ^ {p, s} and the Rayleigh fading SISO ehannel between TX p and TX s 

is denoted as ~ CAA(0,(Tg). Also, the elements of ehannels G {p, s}, are spatially 

eorrelated, henee hmn ~ CM{Qm, R-mn), with Tn,n e {p, s} or 

1 

TTt^Tl G {p, 5]•, (1) 

1 

where Hmn is the symmetric square root of covariance matrix Hmn of channel vector hmn and 
^^mn,w ~ CJ\f {Om 5 Im)- 

Regarding the availability of CSIR, a praetieal seenario is eonsidered, aeeording to whieh 
RX z, i G {p, s}, is aware of direet ehannel ha, while it merely has statistieal knowledge of 
the global uplink ehannel, in the form of eovarianee information. Sinee sueh information is 
slow varying, it ean be available at eaeh of the RXs via a low eapaeity/high lateney feedbaek 
link. Sueh a CSIR formulation is ehosen, beeause standard releases for 4G wireless systems 
require that a given terminal is allowed to report instantaneous CSI to its home BS, however, it 
eannot report sueh information to interfering BSs [[U. 

Sinee SS eonstitutes an essential feature of the investigated hybrid CR system, foeusing on 
seeondary eommunieation, eaeh Medium Aeeess Control (MAC) frame of the SU, that has a 
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CSIR at RX p 


CSIR at RX s 

hpp R,pp Rps R-sp Rss 


^SS R-SS 


RX p Rx s 



Fig. 1. The examined hybrid CR system. 


duration of T time units, consists of a) a SS subframe, the duration of which is r time units, 
followed by b) a Data Transmission (DT) subframe, which lasts for the remaining T — r time 
units. Concerning SS, we choose to apply Energy Detection (ED), since it is characterized by 
low implementation complexity and analytical expressions for the false alarm and detection 
probabilities. Eurthermore, it is assumed that the duration of each MAC frame is such that the 
involved wireless channels remain fixed. 

In what follows, we describe the operation of the studied system, during the SS and DT 
subframes of each MAC frame. 

B. Description of SS phase 

Eocusing on the application of ED for SS, it is assumed that TX s senses the wireless channel 
by sampling the received signal, with a sampling frequency denoted by fg, therefore, SS is based 
on = rfs samples. We define event Hq as the one occurring when the primary system is 
idle, and its complementary event is denoted as "Hi. The received signal at RX s for the n-th. 
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n = 1,..., A^, time instant is expressed, for each of hypotheses "Hq "Hi as 



?7[n], 


ho^/P'pXp[n] + ri[n], if 'Hi, 


if Ho 


( 2 ) 


1 < n < where additive noise, ri[n] is a CSCG, independent, identically distributed (i.i.d.) 
process with ri[n] ~ CAA(0, A^o,o)> Pp denotes the fixed transmit power emitted by TX p and the 
information symbol Xp[n] is selected from a CSCG codebook, i.e., Xp[n] ~ 1) and is inde¬ 

pendent of r][n\. Under these assumptions, when ED is applied, based on a detection threshold, 
denoted by e, e > 0, a closed form expression describing the average (over channel fading) 
probability of false alarm, 'Pf{N, e), as well as an approximation for the average probability of 
detection, Vd{N,e), are derived in Q, [l20ll . These expressions are the following 



( 3 ) 


and 



( 4 ) 


In the following, the DT phase for every MAC frame, is described. 

C. Description of DT phase 

Having described the SS procedure, we now focus on the DT subframe of the secondary MAC 
frame. As explained earlier, the operation of the secondary network during the time intervals 
corresponding to these subframes depends on the obtained SS results. Thus, for the description of 
the received signal during the DT subframe, one needs to discriminate between two SS decision 
cases. 

• Case I: Absence of primary transmissions is detected. We denote this event as 'Hq. Whenever 

such an event occurs, TX s transmits using a power level Pg = Pq. On the other hand, 
RX s employs a receive BE vector w = Wo(hss) G for the detection of the signal 

transmitted by the secondary terminal. 

• Case II: Presence of primary transmission is detected. We denote this event as 'Hi. Whenever 

'Hi occurs, TX s transmits using a power level, Pg = Pi- In addition, RX s employs a BE 
vector, w = Wiifigg) G that is designed taking into account the fact that primary 

activity has been detected. 




In the following analysis, the aehievable instantaneous rate at RX s, regarding the investigated 
system model, is derived. 


D. Rate analysis of the secondary system 

For the determination of the aehievable instantaneous rate of the seeondary system, the signal 
model at the RX side needs to be examined. Using events "Ho and "Hi, that were defined before, 
one ean write the expression for the reeeived signal reaehing RX s, after applying reeeive BF, 
given that event i-ik has oeeurred, as 

Vk = w^hsss/^kXs + CkW^hpss/^pXp + wj^ris, (5) 

k G {0,1}, where Xp,Xs denote the standard eomplex Gaussian signals transmitted by TX p 
and TX s, respeetively, and denotes the noise at RX s. For the latter, it is assumed that 
Us ~ CJ\f {Om, No, s^m)- Variable Ck in © is a diserete Bernoulli distributed RV, that models 
the presenee of primary aetivity during seeondary transmissions. In more detail, RV Ck is equal 
to 1 with probability Pr 

Given ®, the aehievable instantaneous seeondary rate of the examined system model ean be 
expressed as 

TZ = TZq + TZi, ( 6 ) 


where TZk, A: G (0,1}, eorrespond to the rates aehieved in Cases I and II, respeetively. More 
speeifieally, term TZ^ ean be expressed as 


TZk = OLk log2 1 + 




Pk 


N 


+ /5fc logs 1 + 




Pk 


0,s 


Nq. + \w]}h„N P, 


(7) 


where © holds under the assumption that ||iua:|| = 1- Coeffieients ak^jdk in © are defined as 


T’ 7 ~ T 

tto = - Tf: -7^0 (1 “ Pf) 5 /7o = - Tf: - Pi (1 ~ Pd) , and 

T-r T-r 

«1 = — 7^ — PoPf', A = — 7^ — PlPd, 

where Vq = Pr{T-iQ) and Vi = 1 — Pq- Having presented the signal model and the aehievable 
instantaneous rate for the seeondary system, in the following seetion we foeus on the reeeived 
signal model for the primary system. 
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E. Primary system operation mode 

Based on the described operation mode of the secondary network, one can write the expres¬ 
sion describing the received signal at RX p, after applying receive BF, provided that primary 
transmission takes place, as 

Zk = v^hpp^/P^Xp + v^hsps/Kxs + v^rip, if Uk, k e {0,1}, (9) 

where Up stands for the additive CSCG noise received by RX p and v represents the applied 
receive BF vector at RX p, which is assumed to be a vector based on the MRC BF solution, 
thus V = hpp = In our analysis, we assume that rip ~ CM {Om, 

Based on the above described system model, in the following sections we initially investigate 
the QoS, quantified by means of the targeted outage probability, that is achieved for primary 
communication, as well as the achievable average rate of secondary communication. Following 
that, we formulate the problem of optimal SS and reception for the secondary RX, with emphasis 
on the maximization of its achievable average rate, given QoS-based constraints, related to the 
operation of the primary system. 

III. Preliminary Analytical Results 

In this section, closed form approximations describing the outage probability of primary 
communication, as well as the achievable average rate of secondary communication, are derived, 
focusing on a combined CSIR assumption. According to this assumption, the direct channel links 
can be instantaneously available by the RXs, whereas the interference links are merely known by 
their channel covariance matrices. In what follows, an approximation of the outage probability 
of the primary RX is derived in closed form. 

A. Outage probability of primary communication 

An outage event is declared at RX p, when, given that primary transmissions take place, 
the Signal-to-Interference-plus-Noise Ratio (SINR) measured at RX p is below a threshold, 
denoted by 70 . In the following proposition, an analytical approximation for the outage probability 
experienced at RX p, is derived. 

Proposition 1. The outage probability of primary communication, for a hybrid SIMO CR system 
can be approximated as 
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Vout^{l-Vd)J^{Po)+VdJ^{P^), 

where function P{x) is given by 


T(x)^ 




M 


(^pp )'70 

PpXj (R.pp^-j-'Yo^^ 


n“.A,(R„)^n;i., 


I A/g(R,pp) Aj(R,pp) 


and A = E 


^ppR-sp^pp 


can be found in closed form by applying / fZ/l Lemma 3]. 


Proof: The proof is included in Appendix El 


( 10 ) 


( 11 ) 


Having derived an analytical approximation describing the outage probability of primary 
communication, in what follows, we derive a lower bound for the average rate of secondary 
communication, given the instantaneous knowledge of direct channel hgs at RX s. 


B. Achievable average rate of secondary communication 

Given the analysis described in III-Dl the achievable ergodic rate experienced at RX s, condi¬ 
tioned on the knowledge of channel hgs, is given by the following expression 

^hss = ^\hss {^o} + ^\hss ) (12) 

where the occurrence of event Tik, A; = 0,1 is considered for each term E,^hss Hence, this 

leads us to the following analysis: 

1) Occurrence of event TLq: When no primary activity is detected, as a result of SS, RX s 
applies the MRC BF solution such as to maximize the direct signal power, i.e., Wq = hgs = 
ijT^. Also, TX s, in its turn, can transmit with its full available instantaneous power, denoted 
by Ppeak, i-e., Po = Ppeak- In this case, the average secondary rate, conditioned on the knowledge 
of channel hgs at RX s, is given by the lemma that follows. 


Lemma 1. The achievable average rate of secondary communication, conditioned on the in¬ 
stantaneous knowledge of channel h^s at RX s, when event TLq occurs, is characterized by the 
following lower bound 

{"^o} P Co, (13) 


where 


Co = 


(^0 P . 1^0 p 


(14) 













II 


and 


C-ofi — In I 1 + 


Ppeak 11 h s 


N, 


0,s 


, Co,i — In j 1 + 




AT I D ^?s^ps/lss 
1 ^ 0,3 + llh, 112 


(15) 


Proof: The proof is included in Appendix El 


2) Occurrence of event Hi: When event l-Li oeeurs, i.e., when the seeondary system adopts 
the underlay CR approach, the receive BF vector Wi and the transmit power Pi, are system 
parameters which need to be designed. 


Lemma 2. The achievable average rate of secondary communication, conditioned on the in¬ 
stantaneous knowledge of channel hgs at RX s, when event TLi occurs, is given by the following 
expression 


{Pi} 


p I /^i ^ 


(16) 


where 


Cig = In 


1 + 


Pi 

N, 


0,s 



(17) 


and 


Cis = In ( 1 + ■^^\wfhss\'^ 


N, 


0,s 


exp 


'wf (1 m + wA /wf (1 m + Wi 


Pmr.sB'ps'Uli 


El 


Pmr,s^ps'UTl 


— exp 


El 


'UtYP inr^s^ps'UTl J \ p/nr,sBpsUIl / 


( 18 ) 


p 

where pinr,s = is the system Interference-to-Noise Ratio (INR), received at RX s, due to 
primary transmission. 

Proof: The proof is included in Appendix [O ■ 


In the section that follows, an optimization problem is formulated, according to which the 
SS parameters are jointly optimized with the receive BF scheme applied at RX s, with the aim 
of maximizing the conditional (for a given, known instant of channel hgs at RX s) average 
rate of secondary communication, subject to constraints, which are destined to protect primary 
transmissions. 

At this point, it should be noted that, in the remainder of the paper, we will focus on an 
interference-limited CR system, i.e., a system in which interference is the main source of signal 
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dl 


degradation, as compared to noise im With such an assumption, it holds that pinrsRpsiui S> 
Nq^s, hence, assuming that A^o,s = 1, the last term of expression (flSl) asymptotically converges 
to ms eq. (5.1.11)] 


exp 


El 


wfpi^^sRpswi:$>i 


> -7 + ln('m^pinr,sRpsiUi). (19) 


Pinr,sRps^l / V^l Pinr,sRps”^l, 

Consequently, incorporating the high INR assumption, the expectation {"^i} becomes as 
follows 


in' ('D 1 high INR^ 0-1 _ /3l _ 

^\hss {El} -)■ 77777 ^^ 1.0 + 


Quantities Vi q and Vi^i are given by 


ln(2)' 


Eifi — Cifi, 


and 


w\ 


Vi^i = In 


exp 


n Im + wi 


Pinr,slbps^l 

(Im + -^.^ssh^s) 


w\ 


PinrjsR-ps^l 


El 


n Im + ■^hsshf^) wi 


Pinr,s^^ps'^l 


7, 


( 20 ) 


( 21 ) 


( 22 ) 


respectively. 


IV. Problem Formulation 

Having derived a lower bound for the average rate of secondary communication as well as a 
closed form approximation for the outage probability of primary communication, an optimization 
problem can be formulated, the solution of which will lead to a rate-optimal scheme of SS and 
receive BF, with respect to the secondary system, given an outage-based constraint, which aims at 
protecting primary communication from harmful interference. More specifically, the parameters 
that need to be optimized in such a direction, are: a) the SS design parameters, i.e., the sensing 
time, T as well as the ED threshold, e and b) the receive BF vector, Wi, applied at RX s, when 


'Such an assumption is realistic for a CR scenario, as the secondary system can be in the vicinity of the primary, following 
a non-cooperative behavior. 
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event T-ii oeeurs. Henee, the investigated optimization problem ean be mathematieally expressed 
as follows 


maximize 


(PI) 


subjeetto "Pout < "Pout, Vd = Vd, llinill = 1, 

0 < Pi < Ppeak, 0 < r < P, £ > 0, 

where, Pout is the predetermined outage-based eonstraint, and Vd is a targeted average deteetion 
probability for the SS algorithm. 

Solving problem (iPTl) proves to be eomplieated. Thus, we propose to determine the rate-optimal 
BF and SS parameters by solving a simpler optimization problem. The objeetive funetion of the 
new optimization problem is a lower bound of the average rate of RX s, whieh is easier to 
manipulate. This lower bound is: C = Co + As a result, the optimization problem to 

be solved is the following 


maximize C 

€C-^ ^ 1 ,T,£,Pl 

subjeetto Pout < Pout, Vd = Vd, IliUill = 1, 

0 < Pi < Ppeak, 0 < r < P, e: > 0. 


(P2) 


At this stage, we ehoose to divide optimization problem (lP2l) into a number of sub-problems. Fo¬ 
eusing on eaeh sub-problem, one parameter is optimized for given values of the remaining design 
parameters, whieh fulfill the eonstraints. 


V. Solving the Optimization Problem 
A. Determining the transmit power of TX s 

Clearly, the transmit power level. Pi, that maximizes the average rate of the SU, will be 
satisfying the outage eonstraint determined by the primary system, with equality. Henee, one 
needs to solve the following equation 


1 ~ Pd ) P (Ppeak) + PrfP (Pi,root) — P 


out 5 


with respeet to parameter Pi,root- As a result, the following equation is obtained 

Pout - (1 - Pd)-P (Ppeak) 


Pi,root = -p ivo), where yo = 


Vd 


(23) 


(24) 



14 


The inversion of funetion J^(-) leads to a non-closed form expression, thus, a root finding method 
can be applied in terms of solving equation -T'(-Pi,root) — 2/o = 0, with respect to Pi,root > 0. Hence, 
taking into consideration the peak power constraint at TX s, the solution becomes 

P* = min {Pi, root, Ppeak} • (25) 

In the section that follows, an iterative scheme of jointly optimizing the receive BF vector and 
the SS parameters, is thoroughly described. 

B. Jointly optimizing the receive BF vector and the SS parameters 

Having determined the applied transmit power at TX s, P{, which satisfies the outage proba¬ 
bility constraint of problem (lP2l) with equality, the resulting optimization problem that needs to 
be solved is the following 

maximize C 

subject to Vd = Pd, ||'*^i|| = 1, 

0 < r < P, e > 0. 

One can write the objective function of optimization problem (|P^ as follows 

C('u;i, r, £, Pi*) = Co(r, e) + {Pi(mi, r, e, Pj^)} . (26) 

In order to approximate the solution to this problem, we propose to use an iterative procedure 
based on alternating optimization of the SS parameters and the receive BF vector. Following 
such an approach requires solving the following two sub-problems. 

1) Optimizing the SS parameters for a given BF vector: We start with fixing the receive BF 
vector to be an arbitrary unit-norm vector, i.e., Wi = ibi, ||it>i|| = 1. As a consequence, the 
resulting objective function of problem (|P^ is only a function of SS parameters r and e, i.e., 
C = C(iui, r, £, P{). As a result, optimization problem (|P^ becomes 

maximize C(r, e) = ao(r, e)Co,o + /3o(a e)Co.i + ai('r, £)Pi,o + A(p £)Pi,i 

(P4) 

subject to Pd(r, e) = P^, 0 < r < P, £ > 0, 

where terms Pi,o and Pi,i are given by equations (|2T]) and (l22l) . respectively, with Wi = Wi 
and Pi = Pi*. Exploiting the equality constraint for the average detection probability, along 
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with expression dH), one ean express the ED threshold, e, as a funetion of sensing time, r. This 
expression is the following 


e(r) = A^o,o 




(27) 


Substituting (1271) to the objeetive funetion of (IP^ . the following lemma ean be proved, whieh 
is useful for the solution of (|P^ . 


Lemma 3. Function C{t, s^t)) which is obtained after substituting (1271) to the objective function 
of dEH), is a concave function for every r G (0,r]. 

Proof: The proof is included in Appendix [D] ■ 


Sinee the resulting optimization problem is a eonvex problem, any eonvex optimization al¬ 
gorithm ean be applied (i.e., a gradient aseent-based algorithm), with the aim of finding the 
rate-optimal values t* as well as £* (through (1271)). for the given reeeive BE veetor, Wi. 

2) Optimizing the receive BF scheme for fixed SS parameters: The problem of designing 
reeeive BE veetor Wi, sueh as to maximize the objeetive funetion of problem (iPSl) . for given SS 
parameters that satisfy the deteetion probability eonstraint, is equivalently expressed as follows 

maximize {7^i(mi, f, £, )} 

subjeet to VdirJ) = Vd, ||ini|| = 1, 

where f G (0, T] and e > 0. The objeetive funetion of problem (lP5l) is given by (l20l) . with 
Pi = Pj, ai = ai{T,e) = ai and fii = fii{f,e) = fii. Consequently, incorporating the high 
INR assumption, the objective function of the receive BE problem becomes as follows 

\wfR^f{Wi 

(28) 


P* LT 

where Hgff = 1m + j^hsghf^ and Rgff = Pinr,sRps- The lemma that follows assists in solving 
problem (iPSl) with respect to vector Wi, when the objective function is given by expression (l28l) . 

Lemma 4. Considering an interference-limited (high INR) system scenario, optimization problem 
(iPSl) can be approximated by the following problem, the objective of which is a lower bound of 


„ JT? 1 INR^ 


h^ 


In (mf Hgfftui) -f 


A 

ln(2) 


exp 


mfRgffmiy ^\wfRssWi) 
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the objective of problem (lP5l) 



(P6) 


where 



(29) 


and 



(30) 



Proof: The proof is included in Appendix El 


Problem (IP^ . i.e., the problem of maximizing the sum of a quadratic form and a Rayleigh 
quotient over the unit sphere, can be efficiently solved by applying the Trust Region Self 
Consistent Field (TRSCF) algorithm which was introduced and evaluated in lf22l Algorithm 
2 ]. 

An interesting sub-case, which is worth investigating, is the case where ^ 1, Pi —)■ 1. 
When the primary system is (almost) always in transmission mode, then, by focusing on an 
interference-limited system scenario, we obtain an expression reminiscent of [HI eq. (12)] for the 
SIMO interference channel, which is the following 




(31) 


In such a case, the optimal receive BF can be found as shown in the following proposition. 

Proposition 2. For fixed SS parameters, along with a given transmit power level, Pj, which 
satisfies the constraints of (lP2l) . and assuming that Vd —^ f well as that the 

investigated system is interference-limited, the optimal receive BF vector at RX s, in terms of 
maximizing the conditional (with respect to channel hss) average rate of the secondary system, 
is given by the following expression 



where Pj has been obtained in (l25l) . The solution of the latter problem is the eigenvector that 

_2 / p* u \ 

corresponds to the dominant eigenvalue of matrix {pmr,s^ps) ( Im + ’N^^sshfsf 
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Proof: The proof is included in Appendix 0 ■ 

3) Iterative optimization framework: Having solved separately the SS and BF optimization 
problems, we propose to approximate the solution to the rate-optimal joint SS and BF design 
by applying the following iterative algorithm. 


Algorithm 1 Jointly optimizing BF veetor Wi and SS parameters r and e 

1 Initialization in = 0). Fix the reeeive BF seheme sueh that Wi = wf''^ and inerease eounter 
by one. 

2 For the n-th iteration, solve problem (|P^ with Wi = and find values and 

3 Utilizing values r„ and solve problem (IP^ and determine BF veetor 

4 Compute the value of the objeetive r„, e„). 

5 Inerease the eounter by one and if |C„ — C„_i| < where n > 2 and (^ > 0, G M is an 
arbitrary small number, stop, otherwise go to Step 2. 


Remark 1. Sinee for Vi —?• 0, = iumrcj while in the ease where Vi —?• 1, the optimal 

reeeive BF veetor is the DGE of matriees Heff and Reff, i.e., wl = argmaxp||=i = 

iudgej a heuristie ean be exploited in terms of ehoosing veetor wf‘\ For instanee, one ean use 
= iumrc when the primary aetivity profile is low, otherwise veetor = iudge can be 
used. Sueh a heuristie ean be proved useful in terms of redueing the eomplexity of Algorithm 

m 

In what follows, we foeus on the standard interweave and underlay CR systems, and opti¬ 
mization problems, equivalent to (lP2l) . are formulated and then solved. 

VI. Optimizing design Parameters for Standard CR Systems 

The goal of this seetion is to derive rate-optimal system designs for interweave and underlay 
CR systems. In what follows, we start with the interweave (opportunistie) CR system. 
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A. Interweave CR system 


Focusing on the interweave CR system, we assume that TX s transmits with a fixed power 
level, Ppeak and the reeeive BF veetor at RX s is based on the MRC solutioil^. The rate-optimal 
design for the interweave CR system boils down to the following problem 


^ Cllo(Anti S^int) ^ /^O (Ant i £^int) ^ 

maximize Cnt = — 7 - 7 ^. -Co,o H-—Co,i 

-fmt.eint lia(2) 


(P7) 


SUbjeet to 7^out,int — ^out; 0 < Tint ^ £int ^ 0, 
where the objeetive is given by the lower bound in (fT4l) and quantities ao(Ant, Ant) and /3o(Ant, ^int) 
have the form of ao and /3o, (given in ([8])), respeetively, by exploiting the new SS parameters. 
Tint and Sint- The following proposition will be useful for solving problem (IFTI) . 


Proposition 3. The outage probability of primary communication, for a SIMO interweave CR 
system is approximated by the expression that follows 


'Poucint ~ (1 — VdjP{Ppeak) + PdG) 


(33) 


where 


^ AAR 


Q = 


nLv(R 


E 


(Rpp) 


exp 


70 


psnr,p^j (ri-pp) 


nf=i, 


VPJ j = l 


y^fcCR-pp) 


(34) 


and Psnr,p = stands for the system Signal-to-Noise Ratio (SNR) observed at RX p. 
Proof: The proof is included in Appendix 


Having derived a elosed form approximation for the outage probability of primary eommuni- 
eation, one ean express the ED threshold, £int, as a funetion of sensing time. Tint, after substituting 
(l3^ to the outage probability eonstraint of (IFTI) and exploiting the elosed form approximation 
for the average deteetion probability, whieh is given by dH). This expression is the following 

£tnt = <5f^fc= + lV (35) 

Vv Ant/s / 

where, ^int = Q~^ ^ ^ Substituting (1351) to the objee¬ 
tive funetion of (IFTI) . one obtains a single variable objeetive funetion: W(Tint) = -f 


^Regarding the transmit power of TX s for the interweave case, we choose power level Ppeak, because it is assumed that 
Pd 1 . 
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By applying the second derivative criterion, it can be shown that W(Tint) is a concave 
function of its argument, when Tim G (0,r], consequently, an optimal and a corresponding 
(by exploiting ([35])) can be found, by applying a convex optimization algorithm. 

In what follows, the optimal parameter design problem is formulated and solved for an underlay 
CR system. 


B. Underlay CR System 

Concentrating on the corresponding underlay CR system, the optimization problem, equivalent 
to (lP2l) . that has to be solved, is the following 


maximize Cund 

l«ijndGC*fxl,Pund 


(P8) 


subject to ^outjUnd ^ ^outj ||'*^und|| — Ij 

where, lUund represents the applied receive BF vector at RX s and Pund denotes the transmit 
power of TX s. Due to the lack of a SS procedure (r = 0), the conditional average rate of 
secondary communication is given by the following expression 


Cnnd 


' ^ 'T'^und I ' 1 'T^und 
■^1,0 ^1,1 ’ 


(36) 


ln(2) ln(2) 

where quantities P“"q and are given by (|2TI) and (l22l) . respectively, with mi = lUund and 
Pi — Pund- A closed form approximation of the outage probability of primary communication, 
considering an underlay CR system, denoted as Pout,und, is given in the following proposition. 


Proposition 4. The outage probability of primary communication, for a SIMO underlay CR 
system is approximated by the following expression 


V, 


out,und 




und) • 


(37) 


Proof: The outage probability of primary communication is given by 


P, 


out,und 


= Pr 


P \\h 11^ 

p 11 ' ^pp \ I 


< 7o 


Aq^P + Pund\h^ph 

The latter probability has been approximated in Appendix 0 which concludes the proof 


(38) 


Having derived an approximate expression for the outage probability of primary communi¬ 
cation, in closed form, problem (iPSl) can be efficiently solved. More specifically, by following 






20 


the steps of Seetion IV-AI with —)■ 1, the transmit power of TX s can be determined and by 
following the steps of Section IV-B21 with di = Vq and /3i = Pi, one can find the rate-optimal 
receive BF vector at RX s. 

In the following section, the throughput performance of the designed hybrid CR system is 
numerically evaluated and compared to the throughput performance achieved by the designed 
interweave and underlay CR systems. 

VII. Numerical Evaluation 

In this section, the throughput performance of the designed hybrid CR system is evaluated 
and compared to the throughput performance achieved by the equivalent standard interweave and 
underlay CR systems. We use Monte Carlo (MC) simulations with 2500 channel realizations, in 
order to evaluate the performance of the designed CR systems. An interference-limited system 
is assumed, the parameters of which are included in Table HI It should be noted that the values 
of these parameters remain fixed in the remainder of this section, unless otherwise stated. 
The exponential antenna correlation model is adopted, as described in Il23ll . More specifically, 

TABLE I 

Basic simulation parameters 


MAC frame size, T 

100 msec 

Number of receive antennas, M 

4 

SINR threshold, 70 

3 dB 

Sampling frequency, fs 

6 MHz 

Noise variance. No = No,o = No,p = Vo,s 

0 dB 

Antenna correlation factor, p 

0.5 

Power level, Ppeak 

10 dB 

Power level, Pp 

10 dB 

Variance of TX p-TX s channel, erg 

-3dB 

Targeted average detection probability, Vd 

0.975 


considering the (p, g)-th element of the covariance matrix of channel fi-mm m, n G {p, s}, it is 
taken to be p,g = l,...,M, pG[0,l]. 

We start with evaluating the quality of approximating Poutage(P), 7o) = Pr (^ Nq p^Po\h^^hsp\^ ^ 
by value P{Po) for different values of SINR threshold 70, when p = 0.2 as well as when 
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Fig. 2. Quality of approximating probability Poutage(F'o, 7 o) by value J-{Pq) for different values of 70 . 


p = 0.5. As it is evident from Fig. [2l the approximation is satisfactory for the examined range of 
7 o when p = 0.2, but also when p = 0.5. In the latter case, the approximation quality becomes 
higher for relatively high values of 70 . 

In Fig. [3l the average rate of RX s is depicted as a function of the outage probability of 
primary communication, when the primary system is in transmission mode for 30% of the 
time. The throughput performance of the optimized hybrid CR system is plotted together with 
the one achieved by the optimized interweave and underlay CR systems. One can observe that the 
performance of the hybrid CR system overcomes the one achieved by the standard CR systems for 
the whole examined outage probability range. Also, all three curves are monotonically increasing, 
which can be explained by the fact that, as the outage probability constraint becomes looser, the 
secondary system can utilize its available resources primarily with the aim of maximizing its 
spectral efficiency. In addition, the average secondary rate, achieved by the interweave system 
outperforms the one of the underlay system for almost the whole examined outage probability 
interval. This happens, because for low primary activity profiles and for the given quality of the 
SS channel, it is better to sense the existence of spectral “holes” in time, in order to then exploit 
the full potential of the secondary system’s resources (i.e., full transmit power). 

The same performance metric is illustrated in Fig. IH this time for a high activity profile of 
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Fig. 3. Ergodic rate of RX s vs. outage probability of primary communication, Vi = 0.3. 



Fig. 4. Ergodic rate of RX s vs. outage probability of primary communication, Vi — 0.7. 
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the PU, i.e., when it is aetive for 70% of the time. In this ease, the following observations 
can be made: a) The performance of the optimized hybrid system always overcomes the one 
achieved by the optimized standard CR systems, however, the average secondary rates of all 
systems are lower than the ones achieved given a low primary activity profile. This occurs 
because more interference from the primary system is received by RX s, on average, b) The 
underlay CR system now outperforms the interweave one for almost the whole investigated 
outage probability interval. Such behavioral change can be explained by the fact that, as the 
primary system transmits more frequently, it is better for the secondary one to exploit the full 
duration of the MAC frame for DT. 



Fig. 5. Ergodic rate of RX s vs. primary activity profile, Vom = 2 x 10 ^. 

In Fig. [51 the achievable average rate of RX s is depicted for the three investigated systems, as 
a function of the activity profile of the primary system, when the outage probability of primary 
communication is equal to 2%. One can observe that the average throughput of RX s regarding 
the hybrid system, balances between two “extremes” with respect to the activity profile of the 
PU. More specifically, the hybrid CR system behaves similarly to the interweave one, when 
the PU is idle for most of the time, whereas it approaches the throughput performance of the 
underlay system, when the PU is active for most of the time. Also importantly, all three curves 
are decreasing. This occurs because, when the primary system is busy for an increased fraction 
of time, more interference will be received by RX s, on average. 
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Fig. 6. Transmit power of TX s vs. targeted PU outage probability, Vi = 0.3. 

Finally, in Fig. the transmit power levels Pi of the hybrid CR system and Fund of the 
equivalent underlay CR system are shown as a funetion of the targeted outage probability of the 
PU, when the primary system is aetive for 30% of the time. It is elearly shown that the transmit 
power of TX s inereases as the PU beeomes more tolerant to interferenee. It is also observed 
that when the outage probability constraint becomes very loose, i.e., when Pout is about 25%, 
both the hybrid and the underlay systems exploit that flexibility in order to transmit with full 
power. 


VIII. Conclusions 

In this paper, the uplink of an interference-limited, hybrid interweave/underlay CR system 
has been studied. Correlated Rayleigh fading has been assumed for the involved direct and 
interference channels. A realistic CSIR scenario has been examined, according to which each 
RX has mere access to the instantaneous link of its assigned terminal, along with covariance 
information regarding the global uplink channel. Novel closed form approximations, describing 
the outage probability of the primary system have been derived, considering the hybrid CR system 
as well as the standard interweave and underlay CR systems. Exploiting these expressions, 
a new optimization problem that consists in jointly optimizing the SS parameters and the 
applied receive BF scheme, towards maximizing the average secondary rate, subject to an outage 
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probability-based constraint for primary communication, has been formulated and solved for all 
the examined systems. It has been numerically shown that the optimized hybrid CR system 
outperforms the equivalent, optimized interweave and underlay CR systems, in terms of spectral 
efficiency. Interesting extensions can be made regarding the existence of multiple secondary 
terminals. 

Appendix A 

Proof of Proposition □ 

The outage probability of primary communication is given by the following expression 



(39) 


where term Vout,k corresponds to the occurrence of event Hk, k = 0,1. Hence, for the first term 
of (l39l) . one obtains 



(40) 


Let us start by defining the RVs X = h^phgp and Y = |Xp. Assuming that vector hpp is given 

i ~ i 

and that hgp = where ~ CM{Om, Im), one can write X as: X = h^pRiphsp^u,. It 


then holds that, given hpp, A is a complex normal RV with zero mean and variance ax = 


iih“Riii = Mi" . therefore, variance ax is a ratio of quadratic forms. Hence, given hpp. 


Y is an exponential RV with mean value equal to a^, i-O-, it has a Probability Density Function 
(PDF) of the form 



(41) 


and a Cumulative Distribution Function (CDF), given by 



(42) 


Conditioned on the knowledge of hpp, one can write (l40l) as 



( 43 ) 
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Assuming that covariance matrix Rpp has M distinct eigenvalues, RV Z = ||^ppP is distributed 
with PDF given by eq. (14)] 


fziz) 


nf=iA,(R, 


■pp, 


M 

Sn 


exp 


\j (Rpp) 


M 

k=l^k^j 


A/c(R,pp) Aj(Plpp) 


z>0. 


(44) 


As a result, probability Pout,o can be approximated by the following expression 

Voutfl ^ (1 - Vd) [ exp f — fz{z)dz, (45) 

Jo V 7 o^oA PoAy 

where A = The latter expeetation ean be eomputed in elosed form by exploiting ETl 

1 1 

Lemma 3], with matriees A = Rpp and B = RppRgpRpp. Consequently, term Pout.o beeomes 


Pout,o~ (l-Pd)-^(n), (46) 

where funetion T{x) is given in (fTTI) . Following a similar analysis, probability Pout,i is given 
by the following expression 

(47) 


whieh eompletes the proof. 


Appendix B 
Proof of Lemma [I] 

In the oeeurrenee of event I-Lq, the average seeondary rate, eonditioned on the instantaneous 
knowledge of ehannel hgs at RX s, is given by 


{^o} = tto loga (^l + 


P \\h l|2 

^ peak II' ^ss II 


N, 


0,s 


+ /5oiE.|hss ) iog2 ( 1 + 


i"peak||fc..|P 
d^0,s + 1^0 


where Uq = Pp\hf^hps\'^. For the expeetation appearing in the seeond term of 
Jensen’s inequality with respeet to ehannel hps, we obtain 

T^peak II II 


^hss i log2 1 + 


No,s + Uq 


> log 2 1 + 


= loga 1 + 


P 11^ l|2 
peak II II 

No,s + 


P III) l|2 

peak II '^ss II 

AT I D h^gUpshss 

1^0,8 + 11^ 1,2 


(48) 
by applying 


( 49 ) 


whieh eompletes the proof. 















Appendix C 
Proof of Lemma [2] 
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The achievable average rate of secondary communication, conditioned on the instantaneous 
knowledge of channel hgg at RX s and given that event 'Hi has occurred, is given by the following 
expression 


{^ 1 } = «llog2 (^1 + 


Pi\w^h 


N, 


0,s 


+ S log 2 I 1 + 


+ Ul 


(50) 


where Ui = Pp\wfhps\'^. Considering the second term of (l50l) . we get the following expression 


^hss i log2 1 + 


Pl\w^hss\' 
No,s + 


= E 


|h.3 


log2 1 + 


No,, 1^1 

1 + n 


(51) 


where, Yi = Pmv,s\'Wihpsf and BF vector Wi is independent of hps, since no instantaneous 


knowledge of hps is presumed. RV Yi can be written the following way 


) 


(52) 


where hps,^ ~ CA/'( 0 m, Im)- 


It is, thus, easy to confirm that Yi is an exponentially distributed RV, and its PDF is given by: 


frAvi) = 


Pinr,s'*AJ^ rtpsl^l 


exp 


yi 


p[nr,s'^l R’ps'W^l 


. As a result, for the expectation in (ISTl) . one obtains 


^hss { log2 1 + 


No,, 1^1 

1 + n 


= E|I^33 \ log2 ( 1 + Vi + 


(53) 


-E|h„ {log 2 (l +Vi)}. 

For the first term of (l5^ . by exploiting [|2^ eq. (4.337.1)], we obtain the following expression 


E|h„ |log 2 (^1 + Vi + 


= log2 ( 1 + 


N< 


0,s 


ln(2) 


exp 


w. 




El 


? i'^M + ^^hssh^s) 


(54) 


”1^1 Pinr,sRps^l 


Also, exploiting [|2^ eq. (4.337.2)], one can derive the second term of (15^ as follows 


{log2(l + 3"i)} = 


ln( 2 ) 


exp 


Pinr,sRps”I^l 


El 


Pinr,sRps'^l 


(55) 


Substituting expressions 
proof. 


to (l50l) . expression (fT^ is obtained, which completes the 
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Proof of Lemma [3] 
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The resulting single-variable objective function of (lP^ is expressed as 

C(r, s{t)) = aoCo,o + /5oCo,i + aiDi^Q -f 


( 56 ) 


where ai, /Si, i = 0,1 have been defined in ([8]). Taking the derivative of (15^ . with respect to r 
and letting 6 = A^o,o (l + ^^d ^ = Q~^{Vd), we have 


c)C(r, £(r)) _ Vo , r ^ r ]0 ( 

--— —^(6o,o + Lo,i) — —{T>i^o — Lo,o)y ( V^Js 


exp —- 


N< 


^ -lU 


0,0 


N< 


{T — T)Vo{Vifl — Cpfi) 

T^/2^ 

leK)Pi(»i.i - c„,i)^ 


Taking now the derivative of (157]) with respect to r, one obtains 


^0,0 

2^ 


1 + 




iVf 


0,0 


N, 


0,0 


- 1 


0,0 


fs 


2v/^ 

(57) 


d'^C{T,e{T)) _ Vo{A,o-Co,o) 


dV 

{T — yVolVifi — Cpfi) 

tV^ 

{T — r)Po(^i,o — Co,o) 


exp 


1 / 


21^ 




-1 + 


0,0 




N, 


0,0 


N, 


- 1 


0,0 




^ s 




-1 + 




N, 




exp 




N, 


0,0 


-1 + 




- 1 


0,0 


0,0 


<5^ 


N, 


0,0 


^Jl_ 

4tL 

(5 


exp 






-1 + 


0,0 




N, 


0,0 


N, 


0,0 




(58) 


C jj '' 

In the above expression we have — 1 = > 0. Also, Vi^q — Co,o < 0, since it holds 


6 1 Pv'^'^ 

Vo,( 

that < Ppeak and < ||/i 5 s||^. As a result, ^ < 0 and, thus, according to the 

second derivative criterion function, C(t, £(r)), is concave when r G (0,T], which completes 
the proof. 


Appendix E 
Proof of Lemma [4] 

We define variables xi = mfHeffiUi, X 2 = and function /(xi, X 2 ) = /o(a;i) + /i(a: 2 ), 

1^2 xv-eff^l 

where /o(xi) = diln(xi) and /i(x2) = /3i(ln(x2) + exp(x2)Pi(x2)). 

One can easily observe that: (i) function /o(xi) is defined for xi G ^0 — [Amin(Heff), Amax(Heff)] 
and it is concave in the same interval and (ii) function /i(x2) is defined for X2 G = 
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[Amin(Rgg^Heff), Aniax(R'eff Heff)] and it is concavc in the same interval. Sinee the two single¬ 
variable funetions are eoneave within their domains, it ean be eoneluded that 

. 'ixi G Ao, fo{xi) > zo = K 1 X 1 +K 2 , where line zq is defined by points (Amin(Heff), /o(An,i„(Heff))) 
and (Afnax(Heff),/o('^max(Heff))) and 

• yx 2 G All, fi{x 2 ) > zi = fiiX 2 + /i 2 , where line zi is defined by points (Amin (R-eff Heff), 

/l(Amin(R-eff Heff))) and (Amax (R-eff Heff) ,/l (Amax (R-eff ^eff) ) ) • 

As a result, instead of solving optimization problem 

= arg max^/o(m^Heffm) + /i 

for fixed r = f and £ = e, an approximated version of it ean be solved, where the new objeetive 
is a lower bound of the objeetive of problem (IP^ . This completes the proof. 


Appendix F 

Proof of Proposition [2] 

The optimal receive BF vector is obtained by solving the following optimization problem at 
RX s: 


wl= arg max 

11=1 


(P) 


where an approximation of the objective for the investigated regime is given by 
One can rewrite as 


where = 


w^{1m + ) wi 

"IPl Pinr,sRps'^l 


(59) 


(60) 


Focusing on the fact that ^ 1, Pi —)■ 1, function V(-) is defined as 

V(^^J = (In (/2^J + exp(/i^jPi (//^J + 7). 

ln(2) 

By differentiating and using [fT^ eq. (5.1.26)], one can prove that V(/ii„J is an increasing 

function of ■ Consequently, the optimization problem © is equivalent to the Rayleigh - Ritz 
quotient maximization problem 


Wi = 


arg max . 

11 = 1 


(61) 
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By setting the derivative of with respeet to Wi, equal to zero, it ean be found that the 
optimal BF veetor is the one satisfying the equality 


Im “f T7 ^sshgg ] II,u^^P{^j-gRpgW^. 


N< 


0,s 


(62) 


As a result, by inspeeting (l62l) . one ean eonelude that the optimal BF veetor for RX s is the 
dominant generalized eigenveetor (DGE) of matrix pair ( (Im + '^^^sshfg] ,Pmr.sR-ps)• 


Appendix G 

Proof of Proposition [3] 

Foeusing on an interweave CR system, the outage probability of primary eommunieation is 
defined as follows 


^out.int = {I-Vd) Pr 


P \\h IP 

p\\* *'pp \ I 


^0,p -Ppeak I I 


<lo]+V^Pr' 


N, 


< 7 o 


0,p 


A, 


Exploiting the proof steps of Appendix the first term of (1631) . is given by 

P^out,int,l ~ (1 - 'Pd)J'(Ppeak). 

The seeond term of (l63l) . is given by the following expression 

P^out,int,2 PdP^ (Psnr,p|| ^pp|| 7 To) • 

Sinee the PDF of RV Zi = Psnr,p||^pp|P, is known, we obtain 

no 

"Pout,int,2 = Pd fZi{Zl)dZi = VdG, 

Jo 

where Q is given in (l34l) . This eompletes the proof. 


(63) 


(64) 


(65) 


( 66 ) 
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